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ABSTRACT 

We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo 
L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth, and fast- 
dump digital spectrometers. The search targets low Galactic latitudes, \b\ < 5°, in the accessible longitude 
ranges, 32° < I < 77° and 168° < I < 214°. The instrumentation, data processing, initial survey observations, 
sensitivity, and database management are described. Data discussed here were collected over a 100 MHz 
passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 fj,s. Analysis of the 
data with their full time and frequency resolutions is ongoing. Here, we report the results of a preliminary, 
low-resolution analysis for which the data were decimated to speed up the processing. We have detected 29 
previously known pulsars and discovered 1 1 new ones. One of these, PSR J1928+1746, with a period of 69 ms 
and a relatively low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified 
EGRET source 3 EG J 1928+ 1733. Another, PSR J 1906+07, is a non-recycled pulsar in a relativistic binary with 
orbital period of 3.98 hr. In parallel with the periodicity analysis, we also search the data for isolated dispersed 
pulses. This technique has resulted in the discovery of PSR J0628+09, an extremely sporadic radio emitter 
with a spin period of 1.2 s. Simulations we have carried out indicate that ~ 1000 new pulsars will be found in 
our ALFA survey. In addition to providing a large sample for use in population analyses and for probing the 
magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond 
pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams 
from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and 
possibly new classes of transient radio sources. 

Subject headings: pulsars: general — pulsars: individual (PSR J0628+09, PSR J1906+07, PSR J1928+1746) 
— surveys 



1. INTRODUCTION 

Radio pulsars continue to provide unique opportuni- 
ties for testing theories of gra vity and pro b ing states of 
matte r otherwise inaccessible JStairsI 120031 iKramer et alJ 
2004). In large samples, they also allow detailed mod- 
eling of the magnetoi onic components of the int erstel- 
lar medium (e.g., ICordes & Laziol 120021 iH anl 12004] ) and 
the Galactic neutron star population (Lorimer et al.l 1 19931: 
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For these reasons, we have initiated a large-scale pulsar 
survey that aims to discover rare objects especially suitable 
for their physical and astrophysical payoffs. Of particu- 
lar importance are pulsars in short-period relativistic orbits, 
which serve as important tools for testing gravitational the- 
ories in the strong-field regime. Our survey parameters and 
data processing are also designed to find millisecond pulsars 
(MSPs). MSPs with ultrastable spin rates can be used as 
detectors of long-period (> years) gravitational waves (e.g., 
Lommen & Backer! 120011 IWvithe & Loebl 12003b Uenet et all 
2004), while sub-millisecond pulsars (if they exist) probe the 
equation of state of matter at densities significantly higher 
than in atomic nuclei. Long-period pulsars (P > 5 s) and 
pulsars with high magnetic fields are also of interest with re- 
gard to understandin g their conn ection, if any, with magne- 
tars iWoods & T hompson 2005), and improving our under- 
standing of the elusive pulsar radio emission mechanism. Ad- 
ditionally, pulsars with especially large space velocities, as re- 
vealed through subsequent astrometry, will help constrain as- 
pects of the formation of ne utron stars in core-collapse super- 
novae (e.g.. lLai et al]l20"0ll) . Finally, multiwavelength anal- 
yses of particular objects will provide further information on 
how neutron stars interact with the interstellar medium, on 
supernovae-pulsar statistics, and on the relationship between 
high-energy and radio emission from neutron stars. 

The new survey is enabled by several innovations. First 
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is the Arecibo L-band Feed Array 1 (ALFA), a seven-beam 
feed and receiver system designed for large-scale surveys in 
the 1.2-1.5 GHz band. The 1.4 GHz operating frequency of 
ALFA is particularly well suited for pulsar searching of the 
Galactic plane. Lower frequencies suffer the deleterious ef- 
fects of pulse broadening from interstellar scattering, while 
pulsar flux densities typically are much reduced at higher fre- 
quencies. ALFA was constructed at the Australia Telescope 
National Facility (ATNF) and installed in 2004 April at the 
Gregorian focus of the Arecibo telescope. The ALFA fron- 
tend is similar to the 13-beam system used on the Parkes tele- 
scope for surveys of pulsars and Hi. The Parkes multibeam 
(PMB) pulsar survey of the Galactic pla ne (Manchester^yd, 
1 2001 1 . iMorris et alJ [20021 iKramer et aD 120031 iHobbs et al. 
120041: [Faulkner et al 120041) has been extremely prolific, yield- 
ing over 700 new pulsars in the past seven years. Our sur- 
vey will complement the PMB survey in its sky coverage and 
will exploit the much greater sensitivity of the Arecibo tele- 
scope. However, because of the smaller size of the ALFA 
beams compared to the PMB system, many more pointings 
must be done to cover the same area of sky. In addition to 
providing better near-term localization of pulsars on the sky, 
the sensitivity of the telescope greatly decreases the time spent 
per pointing, which results in much better sensitivity to pul- 
sars in compact binary systems without searching a large grid 
of acceleration values to combat binary motion. 

Second, our initial and next-generation spectrometer sys- 
tems have much finer resolution in both time and frequency 
than the spectrometer used with the PMB, increasing the de- 
tection volume of MSPs by an order of magnitude. This 
comes at the cost of an increase in data rate by two orders 
of magnitude (~ 0.3 TB hr" 1 ), requiring substantial computa- 
tional and storage resources for analysis and archival, which 
are now available. For our own use in the early stages of the 
survey as well as for long-term multiwavelength studies, we 
will archive both the raw data and data products from the data 
processing pipeline. 

Large-scale pulsar surveys using the ALFA system have 
been organized through a Pulsar ALFA (PALFA) Consortium 
of which the present authors and others are members. The 
planning and execution of PALFA surveys is a joint effort be- 
tween NAIC and the Consortium to obtain legacy results for 
use by the broader astrophysical community. Similar consor- 
tia have been organized for other Galactic science and for sur- 
veys of extragalactic hydrogen. 

The plan for the rest of this paper is as follows. Following a 
brief description of the ALFA system in §|2] in §|3]we describe 
the technical details and logistics of our survey, including sky 
coverage, data acquisition and processing, sensitivities, and 
archival of raw data and data products. In §|4]we report on 
initial results from preliminary survey observations that have 
so far resulted in the discovery of 1 1 new pulsars. Finally, in 
§ |5J we outline our future plans and expectations for PALFA 
surveys. 

2. THE ALFA SYSTEM 

The ALFA feed horns are arranged in a close packed 
hexagon surrounding a central horn at the Gregorian focus 
of the Arecibo telescope. Orthomode transducers provide 
dual linearly-polarized signals to cooled receivers. The beams 
from the seven feeds are elliptically shaped with equivalent 
circular beam sizes (FWHM) of 3.35 arcmin. Beam centers 

^http : //alf a . naic . edu . | 
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FIG. 1 . — Contours of the telescope gain with the ALFA system averaged 
over its passband (1225-1525 MHz). Contours show the maximum gain for 
a given azimuth and zenith angle from any one of the seven beams and are for 
the average of the two polarization channels. The gain values were measured 
on MJD 53129 at low zenith angles using the extragalactic source 3C286 and 
an assumed source flux density of 14.45 Jy for our measured band. Coutour 
levels are at -1, -2, -3, -6, —9, -12, -15 and -19 dB from the central peak. 
The heavy contour is at the -3 dB level and the next heaviest contour outside 
the 6-beam pattern is for -12 dB. Slices through the centers of the individual 
beam patterns of the seven feeds are also shown at constant azimuth and con- 
stant zenith angles. The equivalent circular beam width (FWHM), averaged 
over all beams, is 3.35 arcmin at 1 .42 GHz i Heiles 2004). 



of the outer six beams fall on an ellipse of size 11.0 arcmin 
x 12.8 arcmin. Efficient coverage of the sky requires that we 
compensate for parallactic rotation of the beam pattern on the 
sky as the telescope azimuth changes. ALFA can be rotated 
relative to the telescope's azimuth arm to accomplish this. We 
note that, because the seven-beam pattern is elliptical, there 
are small offsets of the beams from their ideal positions as the 
feed is rotated. 

Figure ^ shows the measured gain contours for the feed 
systems. The on-axis gain is approximately 10.4 K Jy" 1 at 
low zenith angles for the central beam but is reduced to an 
average ~ 8.2 K Jy" 1 for the other six beams. The sys- 
tem temperature looking out of the Galactic plane ~ 24 K. 
Updated estimates of these system parameters and details 
on their azimuth and zenith-angle dependences, are avail- 
able at http : //alf a . naic . edu/perf ormance Re- 
ceiver signals are transported via optical fiber to intermediate- 
frequency electronics and backend spectrometers in the con- 
trol building. 

Currently, we are using four Wideba nd Are cibo Pulsar 
Processor (WAPP) systems (Dowdetal. 2000) to process 
100 MHz passbands centered on 1.42 GHz for each ALFA 
beam. As used in our survey, the WAPPs compute 256 lags of 
the autocorrelation function for each of two 3-level quantized 
polarization channels; correlation functions for the two chan- 
nels are summed before recording to disk as 2-byte integers at 
64 /is intervals. 
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FIG. 2. — Regions of the Galactic plane to be surveyed with PALFA, taking into account declination limits of Arecibo and restricted to \b\ < 5°. Hatched 
areas indicate regions covered so far in the precursor survey and large filled circles represent newly discovered pulsars. Small dots designate known pulsars. We 
also show bou ndaries of several L-band surveys t hat have been made in or near these regions, including the PMB survey and single-pixel surveys with Parkes 
1 Johnston et al. 1992) and Jodrell Bank 1 Clifton et al. 1992). The Swinburne Intermediate Latitude Survey 1 Edwards et al. 2001) covered the same longitude 
range as the PMB but at latitudes 5° < \b\ < 15°. Arecibo surveys at 0.43 GHz have covered some of our proposed search areas, but to distances much smaller 
than we can reach owing to the limiting effects of interstellar dispersion and scattering. 



Within one year, we anticipate using new spectrometers 
that will process the full 300 MHz bandwidth of the ALFA 
frontend system with 1024 spectral channels. The PALFA 
spectrometers will employ many-bit polyphase filters imple- 
mented on field programmable gate array chips to provide the 
channelization. We expect that mitigation of radio-frequency 
interference (RFI) will be more robust with the new spectrom- 
eter compared to the WAPP's 3-level correlation approach. 
RFI shows a rich diversity in the overall ALFA band. For the 
initial portion of the survey using the WAPPs, we have there- 
fore selected the cleanest 100 MHz portion of the available 
spectrum, centered at a frequency of 1.42 GHz. 



3. ALFA PULSAR SURVEYS 

Previous surveys of the Galactic plane accessible with 
the Arecibo telescope have either b een conducted at the 
lower radio frequency of 0.43 GHz (p ulse & Tayloilll975t 
IStokes et all 119861: iNice. Fruchter & TayloJ 119951) or with 
less sensitive systems at 1 .4 GHz ( Clifton e t alJ 119921 
iManchester et alJ2001l) . Our new survey with the ALFA sys- 
tem promises to probe the pulsar population in the Arecibo 
sky significantly more deeply than the previous surveys. 

3.1. Sky Coverage 

Our long term plan is to conduct comprehensive pulsar sur- 
veys of most of the sky accessible with the Arecibo telescope 
(declinations between -1° and +38°) but with emphasis on 
the Galactic plane (e.g., \b\ < 5°). The survey results dis- 
cussed here are for the inner Galaxy (40° <£< 75°, \b\ < 1°) 
and anti-center (170° <£< 210°, \b\ < 1°) regions. Figure^] 
shows the regions that have been and will be covered close to 
the Galactic plane. The PMB pulsar survey covered the re- 
gion 260° < £ < 50° and \b\ < 5°, i.e., there is some area of 



overlap in the inner Galaxy. Later we will conduct a survey at 
intermediate latitudes up to \b\ ~ 20° to optimize the search 
for relativistic binary systems and MSPs. 

Our strategy for sampling the sky employs two methods for 
maximizing the efficiency and sensitivity of the survey. As 
shown in Figure |5J three adjacent pointings of ALFA are re- 
quired to tile the sky with gain equal to or greater than half 
the maximum gain. Rather than using this dense sampling 
schem e, we have so far adopted a sparse sampling scheme 
iFreirel2003l) that makes only o ne out of three of these point- 
ings. Monte Carlo simulations (|yiemmings & Cordes Il2004t 
Faucher-Giauere & Kaspi] 120041) indicate that sparse sam- 
pling should detect ~ 2/3 of the pulsars in the surveyed re- 
gion. This scheme has the advantage that more solid angle is 
covered per unit time, though much of it at substantially less 
than half the full gain. The sparse sampling approach exploits 
the large sidelobes for the outer six beams, which are ~ 16% 
(-8 dB) of the peak gain centered ~ 5 arcmin from the beam 
axis (see Figure^. The gains of these sidelobes are approxi- 
mately 0.7 and 1 .6 times the on-axis gains for the Green Bank 
Telescope and the Parkes 64-m telescope, respectively, and 
thus provide significant sensitivity. 

Later on, we will make the two additional passes needed to 
achieve dense coverage. Despite the smaller numbers of new 
pulsars expected in these subsequent passes, they will yield 
more pulsar discoveries than if we were to extend the sparse 
coverage to higher Galactic latitudes, where the pulsar density 
decreases rapidly. 

We are also exploring a multiple-pass strategy, where a 
given sky position is observed two or more times. This ap- 
proach is motivated by the fact that pulsar flux densities are 
highly variable (over and above the fundamental pulsation), 
sometimes by more than an order of magnitude, due to a num- 
ber of intrinsic and extrinsic causes, including nulling and 



4 



Cordes et al. 



(3 

B 



10 


- 


1 1 1 






1 




' 1 1 - 




- 




1 






■ 




- 




2 




2 




- 


5 


















- 


( 1 




1 




1 ] 


- 
■ 








3 




3 




- 







( 2 




2 




C 1 






- 




\ 




I 




- 






( 3 




3 




3 ) 


- 


p, 
















— o 


_ 




2 




2 




_ 








3 




3 






10 




III, 










, , 1 , - 




-10 


-5 









5 


10 



Azimuth Angle (arcmin) 

FIG. 3. — ALFA beam locations on the plane of the sky showing the half- 
power beam widths for three pointings derived from the data presented in 
FigurefTI The pointings are labeled 1, 2 and 3 in a dense sampling grid that 
covers nearly all of the solid angle with at least half the gain of the relevant 
beam. Sparse sampling consists of making only one of these pointings. Note 
the ellipticity of the beams and their pattern for a given pointing. 



mode changes common to many pulsars, eclipses, and inter- 
stellar scintillation (both long-time-scale refractive scintilla- 
tion and fast, diffractive scintillation). RFI is also episodic. 
Our simulations suggest that pulsars can be missed in single- 
pass strategies but that any improvement from multiple pass 
approaches depends on the details and prevalence of flux 
modulations. We are in the process of comparing different 
strategies while also using simulations to fully optimize our 
usage of ALFA. 

3.2. Data Analysis 

To maximize the pulsar yield and overall science return 
from the PALFA survey, we are processing the data twice. 
During the observations, incoming data are transferred to the 
Arecibo Signal Processor 2 , a computer cluster that processes 
the data in quasi-realtime after reducing the time and fre- 
quency resolution to increase throughput. This "quicklook" 
analysis, described below, is primarily sensitive to pulsars 
with P > 30 ms, which are expected to make up the bulk of all 
discoveries. We are currently developing an offline data anal- 
ysis scheme which retains the full resolution of the data and 
will be sensitive to MSPs and pulsars in short period binary 
systems as well as to pulsars with large values of dispersion 
measure (DM). In addition to using a number of different pul- 
sar search codes and algorithms, this latter analysis pipeline 
will also take advantage of the multiple beam data acquisition 
for RFI excision. Analysis of the raw data will be done on 
several computer clusters at the home institutions of members 
of the PALFA Consortium. Further details will be published 
elsewhere. 

The quicklo ok pipeline use s freely available pulsar data 
analysis tools (Lorimer 2001) to unpack and transform the 
correlation functions from the WAPPs to spectra with 256 



channels every 64 (is. The data are decimated in frequency 
and time by factors of 8 and 16, respectively, to allow quasi- 
realtime processing. The resulting data sets with 32 frequency 
channels and 1024 (is time resolution are then corrected for 
the effects of interstellar dispersion by appropriately delay- 
ing low frequency channels relative to the highest one. This 
process is carried out for 96 trial values of DM in the range 
0-980 pc cm -3 . The step size in DM is approximately 2, 4, 8, 
16 and 32 pc cm" 3 with changes in step size at approximately 
62, 124, 253, and 506 pc cm" 3 . Two different searches are 
then carried out on the resulting dedispersed time series, one 
for periodic signals and a second for isolated pulses. 

The analysis for periodic, dispersed pulses follows most 
standard pulsar search schemes (see, e.g., Lorimer & Kramer 
2005), and the software used for t his analysis is based on 
code developed for an earlier survey ( Lorime r et all200(J) . In 
essence, the procedure is to look for harmonically related sig- 
nals in the amplitude spectrum (the magnitude of the Fourier 
transform) of each dedispersed time series. To increase sen- 
sitivity to signals with narrow duty cycles, which have many 
harmonics in the Fourier domain, the amplitude spectra are in- 
coherently summed so that up to 2, 4, 8 and 16 harmonics are 
combined. A list of candidates with signal-to-noise (S/N) ra- 
tios above 8 is then formed and the data are folded in the time 
domain to produce a set of diagnostic plots of the form shown 
in Figure 0] To simplify the viewing of the search output, a 
web-based browsing system was developed for examination 
of candidate signals during an observing session. The most 
promising pulsar candidates are filed for future observation 
and follow up. 

In parallel with the periodicity analysis, we also search for 
isolated pulses in t he 96 dedispersed tim e serie s based on 
code developed bv iCordes & McLaughlin (2003). In brief, 
threshold tests are made on each dedispersed time series after 
smoothing it by different amounts to approximate matched- 
filter detection of pulses with different widths. In addition, 
we consider events defined by clusters of above-threshol d 
samples in a "friends-of-friends" algorithm (Frederic 1995). 
For each pointing, diagnostic plots similar to those shown 
in Figure [5] are generated for inspection within the candidate 
browsing system. This example shows data for the 1.2 s pul- 
sar PSR J0628+09, which we discovered in our s ingle -pulse 
search through its bright individual pulses (see jj4.4l i. The 
three panels for each of the seven ALFA beams show, from 
left to right: events above a threshold S/N > 5 vs. time and 
DM channel; a scattering plot of DM channel vs. S/N; and a 
histogram of S/N. Events appear with largest S/N in the DM 
channel that best matches the pulsar's DM; they also appear 
in neighboring DM chan nels with lower S/N that depends on 
the pulse width ( ICordes & McLaughlinl2003l) . 

3.3. Database Management and Archiving 

PALFA survey results are archived in a MySQL 3 database 
system, which stores sky coverage and data quality informa- 
tion along with results from the data processing. This system 
is also designed to record the results from several different 
search codes that implement the processing system described 
above to allow comparison and optimization. The MySQL 
database includes an observational table with fields that char- 
acterize specific telescope pointings, the resulting raw data 
files, and ancillary information about the observations. An- 
other table reports the results of the preliminary, quicklook 



2 |http : //astron . berkeley ■ edu/$\sim$dbacker/asp ■ html| 3 http://www.mysql.com 
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FIG. 4. — Sample periodicity search output from the quicklook analysis 
showing the discovery of the 69-ms pulsar PSR J 1928+ 1746. Top left: a 
coarse version of the dedispersed time series used to assess basic data qual- 
ity. Top right: S/N as a function of trial DM. The points show detections in 
the periodicity search while the curve is the theoretically expected response 
given the system parameters and pulse width. Middle panels: gray scales 
showing pulse intensity as a function of sky frequency (quantized into fre- 
quency bands) and observing time (quantized into subintegrations). Bottom: 
folded pulse profile obtained by integrating over the whole pass band. 



data analysis, with fields that describe candidate signals and 
whether they correspond to known pulsars or not. There is an- 
other table that tracks the content and location of portable disk 
drives used to transport data from the observatory to process- 
ing sites. Additional tables report results from the data pro- 
cessing that uses the full data resolution, a list of refined pulsar 
candidates, and the status of confirming and other follow-up 
observations. 

A data archive is under development at the Cornell Theory 
Center and will include the original raw data as well as anal- 
ysis products and database mining tools, accessible through a 
data gateway. 4 

3.4. Search Sensitivity 

Data taken in our preliminary survey dwell on particular sky 
positions for 134 s for the inner Galaxy and 67 s for anticen- 
ter directions. The minimum detectable flux density S m i n for 
PALFA with these parameters is a factor 1.6 smaller than for 

-1/2 

the PMB survey, implying a maximum distance D max oc S mia 
about 1.3 times larger for long-period pulsars. The sampled 
volume on axis is accordingly about a factor of two larger for 
long-period pulsars. In our full-resolution analysis, the vol- 
ume increase is even larger for small periods owing to the 
smaller PALFA channel widths and the shorter sample inter- 
val. For P < 10 ms, the searched volume increase can be a 
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FIG. 5. — Single-pulse search output for a follow-up observation of the 
1.2-s pulsar J0628+09. Each row shows data collected by one of the seven 
beams during the pointing. From left to right the plots show: scatter plot of 
events with S/N > 5 vs. time and dispersion-measure (DM) channel; scat- 
ter plot of DM channel and S/N for events; and the number of pulses ver- 
sus S/N. Individual pulses from J0628+09 are clearly seen only in beam 5. 
The distribution of events vs. DM depends on the pulse shape and width 
ICordes & McLaughlin 2003). 



factor of 10 or more. 
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FIG. 6. — Theoretical minimum detectable flux density (S m j n ) vs. P for 
different values of DM. Short dashed lines: for the coarse-resolution PALFA 
data analyzed with the quicklook software that led to the discoveries reported 
in this paper. Long dashed lines: for the Parkes multibeam survey, which used 
96 channels across 288 MHz and 250 /is sampling for scan durations of 2100 
s. Solid lines: for full-resolution PALFA data. For each set of curves, DM 
values from the lowest to the highest curve are 1, 200, 500 and 1000 pc cirT 3 . 
The breakpoint at P ~ 10 ms for the solid curves occurs because we assume 
that the intrinsic pulse duty cycle scales as P~'/ 2 with a maximum of 0.3, 
which occurs at this period. Above 10 ms, the number of harmonics con- 
tributing to detections increases from 1 to 16 (the maximum searched) as the 
duty cycle gets smaller. A threshold of 10cr is used. 
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Figure [5] shows idealized plots of S m [ n versus P for four 
values of DM (from 1 to 10 3 pc cm" 3 ) using post-detection 
dedispersion followed by a standard Fourier analysis with har- 
monic summing. The values of S m i n include the effects of 
radiometer noise and pulse smearing from instrumental ef- 
fects combined with dispersion and scattering in the inter- 
stellar medium using t he "NE2001" electron density model 
dCordes & Laziol2002l) . Scattering has been calculated in the 
model for the particular direction b = 0° and I = 40°. Direc- 
tions at higher latitudes will show less scattering and better 
sensitivity for large values of DM. The results are also based 
on the assumption that pulse amplitudes are constant over the 
observation time that spans many pulse periods, which is ob- 
viously an idealization. Our calculations for the PMB sur- 
vey do not include high-pass filtering in both hardware and 
softwa re that degrades the s ensitivity to long-period pulsars 
(see iMancheste ret al.ll2.001l) . No high-pass filtering is done 
in our analyis, either in hardware or in software. We empha- 
size that the curves in Figure[6]should be interpreted as lower 
bounds on the true values of S m \ D because real-world effects 
such as RFI and receiver gain variations will raise the effec- 
tive threshold of the survey. Our quicklook analysis described 
above, which analyzes data after decimation in time and fre- 
quency, has detection curves about 60% more sensitive than 
those for the PMB survey except for P < 10 ms, for which the 
large sampling time of the quicklook analysis significantly de- 
grades the sensitivity. 



4. INITIAL RESULTS 

For our preliminary survey observations carried out be- 
tween 2004 August and 2004 October, we have used 17.1 
hours of telescope time for 919 pointings in the Galactic anti- 
center and 32.2 hours for 865 pointings in the inner Galaxy, 
covering 15.8 deg 2 and 14.8 deg 2 in each region respectively. 
These numbers and their graphical presentation in Figure 
were obtained using the MySQL database described earlier. 
Table [2 lists the 1 1 new pulsars that we have found so far. 
Table |2] lists the detection statistics of 29 previously known 
pulsars also seen in the quicklook analysis pipeline. Not in- 
cluded here is a detection of the 1.55-ms pulsar B1937+21, 
which was undetected due to the coarse time resolution of the 
quicklook pipeline. The pulsar was, however, easily detected 
when the raw data were folded at their full resolution. The 
high time-resolution data pipeline mentioned above will al- 
low detection of any MSPs missed in the quicklook analysis. 

4. 1 . General Remarks 

Four of the pulsars discovered in the inner Galaxy 
(J2009+33, J2010+32, J2011+33 and J2018+34) are in the 
northernmost region detectable from Arecibo; all of these 
have DMs between 220 and 350 pc cm" 3 . Some, if not all, of 
these objects are possibly associated with the Cygnus region, 
where a large number of supernovae will have produced many 
relatively young pulsars (e.g., Vlemmings et al. 2004 and ref- 
erences therein). Another four objects (J1901+06, J1904+07, 
J 1905+09 and J 1906+07) are in the southernmost region vis- 
ible from Arecibo, where the density of pulsars is known to 
be higher. This region was previously covered by the PMB 
survey, suggesting that our precursor survey indeed already 
surpasses the depth of the PMB survey when a conventional 
pulsar search analysis is done. Subsequent to our discovery, 
J 1906+07 was identified in the acceleration search output of 
the PMB data (Lorimer et al. 2005, in preparation). 

All previously known pulsars were detected in our point- 
ings if they were within one beam radius of one of the ALFA 
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7. — Regions of the Galactic plane surveyed with PALFA to date 
showing the Galactic anti-center region (top) and inner Galactic plane (be- 
low). Dots denote the pointing centers of each 7-beam cluster, filled circles 
show newly discovered pulsars while the open circles designate previously 
known pulsars. 



beams. In addition, we detect some strong pulsars several 
beam radii from the nearest beam center. A coarse analysis 
suggests that our detection rate is consistent with what we ex- 
pect from the sparse-sampling strategy discussed earlier based 
on simulations. A detailed analysis will be done as we con- 
tinue the survey. 

The single-pulse search analysis is notably successful in de- 
tecting 6 out of 1 1 of the new pulsars and 21 out of 29 of the 
known pulsars. These statistics are consistent with the fact 
that the known pulsars tend to be stronger than our new de- 
tections. The single-pulse analysis is valuable both for cor- 
roborating candidate detections from the periodicity analysis 
and, as we have shown in the case of J0628+09, for identify- 
ing pulsars that are missed in the periodicity search owing to 
the intermittency of their pulses. 

4.2. PSR J 1928+ 1746 

The first pulsar discovered in the ALFA survey, 
PSR J 1928+ 1746, also has the shortest period among 
the pulsars discovered so far, P = 68.7 ms. While time 
for follow-up observations on this and the other pul- 
sars discovered has so far been limited, we have made some 
multi-frequency and timing observations of PSR J 1928+ 1746. 
Using the TEMPO 5 software package to analyze 83 arrival 
times from PSR J 1928+ 1746 spanning a 257-day baseline, 
we obtain results presented in Table [3] The timing model 
implies that PSR J 1928+ 1746 is a young isolated pulsar 
with a characteristic age r c = P/2P = 82kyr, a surface 

magnetic field strength B = 3.2 x 10 19 VPP G = 9.6 x 10 11 G 
(assuming a dipolar field) and spin-down energy loss rate 
E = mn = 1.6 x 10 36 / 45 erg s" 1 (where fi = 2ttP~ 1 and / 45 is 
the moment of inertia in units of 10 45 g cm 2 ). 



http://pulsar.princeton.edu/tempo 
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side is accurate to approximately 20%. 



Multifrequency observations from 1.1 to 9 GHz shown in 
Figure [8] suggest that the radio spectrum is nearly fiat, S„ oc 
j,+o.2±o.3 rpjjg quoted error reflects empirical departures from 
the fit and thus includes any systematic calibration errors or 
random errors from scintillations. Estimates of the flux den- 
sities are coarse because we have simply scaled the signal-to- 
noise ratios of the average pulse amplitudes and used typical 
values for the gain and system temperature. The flux densi- 
ties at the higher two frequencies are likely to be influenced by 
modulations from interstellar scintillation (based on DM and 
the likely distance). High-frequency surveys are naturally bi- 
ased toward the discovery of objects with flatter spectra than 
surveys at lower frequen cies. Also, you ng pulsars appear to 
have flatter spectra tLorimer et alJI 1995h . so high-frequency 
surveys of the Galactic plane will be less biased against them. 
PSR J 1928+ 1746 appears to be a prototype fiat-spectrum ob- 
ject of which we can expect to find more in our survey. 

As shown in Figure |9j PSR J 1928+ 1746 lies well within 
the localization map for the unidentified EGRET source 
3EG J1928 +1733. The EG RET source shows significant vari- 
ability (iTorres et alJ2 001) that is indicative of a blazar, but has 
a photon index, T = 2.23 ±0.32, not inconsistent with those of 
known pulsars. If PSR J1928+1746 is the radio pulsar coun- 
terpart to 3EG J1928+1733, the implied efficiency for conver- 
sion of spin-down energy into gamma-rays is ?y 7 = L 7 /£ = 
22%f2 7 (<f/6kpc) 2 , where Q 7 is the solid angle (sr) swept out 
by the pulsar's beam and a photon index of -2 is assumed 
for the gamma-ray spectrum. While the nominal efficiency is 
higher than that of any of the confirmed gamma-ray pulsars 
jThompson et alJll999h . we note that the above calculation is 
strongly dependent on the uncertain beaming fraction and on 
the DM-derived distance to PSR J 1928+ 1746 of 6 kpc. Also, 
the flux measur ement used to calcu late the efficiency from 
the 3EG catalog ( Hartma n et al.l 19991) is the largest (and most 
significant) value, so the implied efficiency should be viewed 
as an upper bound. Two other young pulsars re cently discov- 
ered within EGRET erro r boxes, J202 1+365 1 ( Rober tset al.l 
and J2229+6114 (Hal pern et all 120011) have similarly 
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FIG. 9. — EGRET localization test statistic of the high-energy gamma-ray 
source 3EG J1928+1733 IHartman et alJ1999l) . Contours delimit probabil- 
ity regions, from innermost to outermost, of 50%, 68%, 95%, and 99% for 
the location of the gamma-ray emission. In addition to PSR J1928+1746 
which lies close to the center, we also show the next nearest pulsar, 
J1930+1852, previously advanced iCamilo et al. 2002) as a possible coun- 
terpart to 3EG J1928+1733 even though it lies outside the 99% contour. 



high inferred efficiencies. These pulsars will be excellent fu- 
ture targets for the Gamma-ray Large Area Space Telescope 
(GLAST). 

4.3. PSRJ1906+07 

PSR J 1906+07 was initially attributed a spin period half of 
its actual value owing to the presence of its interpulse. Very 
recently, this object has been found to be a binary pulsar by 
comparing the parameters in Table[Jwith entries in the PMB 
database and by making new observations at Jodrell Bank. 
The orbital period is 3.98 hr, the eccentricity is 0.085 and 
its projected orbital semi-major axis is a\ sin / = 1.42 light s, 
yielding a mass function of 0. 1 1 M Q . Discovery of the binary 
nature of this source and a discussion of its properties may be 
found in Lorimer et al. (2005, in preparation). 

The discovery of PSR J 1906+07 underscores the power and 
great potential of PALFA surveys for finding binary pulsars. 
The large sensitivity of the Arecibo telescope allows us to use 
integration times short enough that many binaries will be de- 
tected without recourse to acceleration searches, as has been 
necessary with the PMB survey. Acceleration searches lead to 
a greater number of statistical trials which, in turn, require a 
higher detection threshold to minimize the number of "false- 
alarm" detections. 

4.4. PSRJ0628+09 

Many pulsars show large modulations of the pulsed flux 
density. In some cases, it is easier to d etect single 
pulses than the periodicity in a Fourier analysis iNicelll99% 
iMcLaughlin & Cordesl2003T) . In our quicklook analysis, most 
of the known pulsars and three of the pulsars discovered in 
the periodicity search (J0540+32, J1904+07 and J2011+33) 



8 



Cordes et al. 



appear strongly in the single pulse search and less so for three 
other pulsars. Some single pulses of PSR J1904+07 are de- 
tected with S/N ~ 36, more than double that found by the 
periodicity search. PSR J0628+09 was discovered only by de- 
tection of its very sporadic single pulses, some of which have 
peak S/N~ 100. The discovery data set had only three large 
pulses in a 67 s scan, too few to allow the pulsar's detection in 
the periodicity search. A periodicity of 2.48 s was determined 
from the arrival times of those pulses. Subsequent observa- 
tions with a greater number of strong pulses and an above- 
threshold detection in the periodicity analysis have allowed 
us to determine the true period of 1.24 s. 

The discovery of PSR J0628+09 clearly demonstrates the 
importance of single pulse searches. As shown in Figure |5j 
these searches are enabled by the simultaneous measurements 
in multiple beams, which allow discrimination between RFI 
and celestial events. Extrapolating from the present sample to 
the whole survey, we can expect to find a significant number 
of pulsars through their single-pulses and not through their pe- 
riodicity. The analysis also may detect radio transients from 
non-pulsar objects, a plausible outcome given the recent dis- 
covery of a transient radio source in the direction of the Galac- 
tic center (Hv man et alJ20 05) and a number of other transient 
radio sources found in a single-pulse analysis of the PMB sur- 
vey (McLaughlin et al. 2005, in preparation). 

5 . FUTURE PLANS AND EXPECTATIONS 

We have described the initial stages of a large-scale sur- 
vey for pulsars using ALFA, the seven-beam system at the 
Arecibo Observatory that operates at 1.4 GHz. Our discov- 
ery of 1 1 pulsars from precursor observations — using a pre- 
liminary data acquisition system that sampled only 1/3 of the 
available bandwidth followed by a quick-look analysis — is 
extremely encouraging. A new spectrometer that uses the full 
bandwidth will become available within the next year. The 
full data processing pipeline, now under development, will 
have excellent sensitivity to MSPs and is expected to yield 
further pulsar discoveries in our existing data. This pipeline 
will include masking of RFI in the frequency-time plane prior 
to dedispersion, a new matched-filtering search algorithm for 
events that have a broader range of frequency-time signatures 
than those encountered for pulsars, and compensation for ac- 
celeration in binary systems. 

In the near future, we expect to begin regular timing pro- 
grams on several telescopes to obtain precise determinations 
of the spin and astrometric parameters of these pulsars and 
others that will be discovered. 

The full survey will take more than five years, depending 
largely on allocation of telescope time. Numerical models of 
the pulsar population, calibrated by results from the PMB sur- 
vey and incorporating measured characteristics of the ALFA 



system, suggest that as many as 1000 new pulsars will be dis- 
covered. 

The raw data from the search as well as the data products 
from the search analysis will be archived and made avail- 
able to the broader community via a web-based portal. Ini- 
tially, the database system will enable our own mining of the 
data for new pulsars and perhaps other astrophysical signals. 
Later, we expect the system to provide opportunities for multi- 
wavelength searches, such as identification of radio counter- 
parts to X-ray sources or to candidate gamma-ray pulsars seen 
in data from GLAST. 
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TABLE 1 

Pulsars Discovered in the PALFA Precursor Survey 



PSR RA(J2000) DEC(J2000) P DM (S/N) SP? Comments 

(ms) (pc cm -3 ) 



J0540+32 


05 h 40 m 38 s 


+32°02'19" 


524 


120 


36 


Y 


Strong, sporadic single pulses 


J0628+09 


06 h 28 m 49 s 


+09° 09' 59" 


1241 


88 




Y 


Discovered as S/N=40 single pulses 


J1901+06 


19 h 01 m 36 s 


+06° 09' 36" 


832 


162 


14 


Y 




J 1904+07 


19 h 04 m 09 s 


+07°39'41" 


209 


275 


15 


Y 


Strong, sporadic 


J 1905+09 


19 h 05 m 16 s 


+09°01'22" 


218 


452 


14 


N 




J 1906+07 


19 h 06 m 5P 


+07°49'01" 


144 


217 


11 


N 


Interpulse; original detection 

at 72 ms; binary with P or t, = 3.98 hr 


J1928+1746 


19 h 28 m 43 s 


+17°46'23" 


69 


174 


19 


N 


First ALFA pulsar; flat spectrum 


J2009+33 


20 h 09 m 39 s 


+33°25'58" 


1438 


254 


13 


N 


Sporadic 


J2010+32 


20 h 10 m 21 s 


+32°30'22" 


1442 


350 


23 


N 




J2011+33 


20 h ll m 47 s 


+33°21'49" 


932 


300 


30 


Y 


Sporadic 


J2018+34 


20 h 18 m 54 s 


+34°32'44" 


387 


226 


24 


Y 





NOTE. — RA and DEC are the right ascension and declination for the center of the beam where the pulsar was found. Typical half-width uncertainty 
in pulsar position is one beam radius (about 1.6 arcminutes) in both coordinates, except for PSR J1928+1746. DM is the DM value at which the search 
algorithm identified the pulsar with maximum S/N. (S/N) is the signal-to-noise ratio of the averaged pulse shape. SP denotes whether individual pulses from 
this pulsar were detected in the single-pulse search. 
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TABLE 2 

Previously Known Pulsars Detected in the PALFA Precursor Survey 



PSR 


P 

(ms) 


DM 

(pc cm" 




DM 

3 ) 


Si 400 
(mjy) 


(S/N) 


A0 
(') 


SP? 


J0631+1036 


287 


125 


148 


0.8 


76 


6.7 


Y 


J 1 OJJTUJU / 


845 


403 


410 


97 


40 


3.2 


Y 


R1 8SQ4-07 


644 


253 


282 


0.9 


42 


2.3 


Y 


R 100^-4-07 


648 


245 


226 


1 ,o 


11)1 


u.o 


Y 


D 1 7U4TUU 


267 


473 


508 


1 7 


82 


2.4 


Y 


J 1904+0800 


263 


439 


424 


0.36 


17 


2.0 


N 


J 1 71/JTW 1 U 


990 


258 


^83 


0.5 


41 


1 .8 


Y 


TlQ0fv4-0QI? 

J 1 7l"JTu7 1 i. 


775 


265 


240 


0.32 


12 


5.4 


Y 


T 1 907+0740 

J 1 ;/W / i V7 / 


557 


332 


353 


0.41 


20 


2.3 


Y 


T 1 Q07+0Q 1 8 


226 


358 


353 


0.29 


18 


4.7 


N 


B1907+10 


284 


150 


198 


1.9 


57 


1.9 


Y 


J 1 7UOTU / J t T 


212 


1 [ 


46 


0.54 


1 3 


1 1 


Y 


J 17UOTU7V// 


336 


468 


452 


0.22 


48 


1 .7 


N 


Tig 10+0714 

J 171 \7 l \7 / 1^ 


2712 


124 


106 


0.36 


14 


1.8 


Y 


R1Q13+1 


404 


242 


240 


1.3 


30 


4.4 


Y 


J1913+1000 


837 


422 


452 


0.53 


26 


1.7 


Y 


J 171 JTlUl 1 




179 


1 70 


0.50 


10 


2.7 


N 


B 19 14+13 


282 


237 


219 


1 2 


150 


1 .8 


Y 


R 1 Q 1 *>+1 ^ 

D171JT1J 


195 


95 


103 


1 9 


74 


2.3 


Y 


D171 UT 1 *T 


1181 


27 


28 


1 .0 


1 8 


3.0 


Y 


B1919+14 


618 


92 


74 


0.7 


41 


0.5 


Y 


B1921+17 


547 


143 


177 




13 


3.0 


N 


B1925+188 


298 


99 


166 




19 


1.9 


N 


B 1929+20 


268 


211 


205 


1.2 


19 


3.0 


N 


B 1952+29 


427 


8 


18 


8.0 


88 


3.7 


Y 


J 1957+2831 


308 


139 


163 


1.0 


54 


1.6 


Y 


B2000+32 


697 


142 


184 


1.2 


43 


2.2 


Y 


J2002+30 


422 


196 


184 




24 


1.2 


N 


B2002+31 


2111 


235 


197 


1.8 


88 


3.3 


Y 



NOTE. — Pulsar parameters P, DM and S1400 are from the ATNF pulsar database 1 Manchester et al. 2005 ). 
DM is the DM value at which the search algorithm identified the pulsar. (S/N) = signal-to-noise ratio of the 
averaged pulse shape. A8 is the angular distance from the nearest beam centroid in which the pulsar was 
detected. SP denotes whether individual pulses from this pulsar were detected in the single-pulse search. 



TABLE 3 

Observed and derived parameters for PSR J1928+1746 



Parameter Value 



Right ascension (J2000) 19 h 28 m 42?48(4) 

Declination (J2000) 17°46'27(1)" 

Spin period, P (ms) 68.728784754(1) 

Period derivative, P 1 .3209(5) X 10"" 

Epoch (MJD) 53448.0 

Dispersion measure, DM (pc cnT 3 ) 176.9(4) 

Flux density at 1400 MHz, S 1400 (mjy) 0.25 



Surface magnetic field, B (Gauss) 9.6 X 10" 

Characteristic age, r c (kyr) 82 

Spin-down luminosity, E (ergs s~') 1.6 X 10 36 

DM Distance (NE2001), D (kpc) -6 

Radio luminosity at 1400 MHz, S1400O 2 (mJykpc 2 ) ~ 9 



NOTE. — Since the timing data collected so far span only 257 days, the phase-connected timing solution should be viewed as preliminary. The figures in 
parentheses give the uncertainties in the least-significant digits quoted. To be conservative, these are calculated by multiplyin g the nominal 1 a TEM PO standard 
deviations by an ad-hoc factor of 10. The DM distance is calculated using the NE2001 electron density model for the Galaxy 1 Cordes & Lazio 2002 ). 



